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Abstract 

Ehlers and Kundt have provided an approximate procedure to 
demonstrate that gravitational waves impart momentum to test par- 
ticles. This was extended to cylindrical gravitational waves by Weber 
and Wheeler. Here a general, exact, formula for the momentum im- 
parted to test particles in arbitrary spacetimes is presented. 

1 Introduction 

There has been debate whether gravitational waves really exist [1,2]. To 
demonstrate the "reality" of gravitational waves, Ehlers and Kundt [1] con- 
sidered a sphere of test particles in the path of a plane-fronted gravitational 
wave and showed that a constant momentum was imparted to the test par- 
ticles. A similar discussion was given by Weber and Wheeler [2,3] for cylin- 
drical gravitational waves. An operational procedure embodying the same 
principle has been used [4] to express the consequences of relativity in terms 
of the Newtonian concept of gravitational force. The pseudo-Newtonian (^N) 
"gravitational force" is defined as the vector whose intrinsic derivative along 
the separation vector is the maximum tidal force, which is given by the accel- 
eration vector for a preferred class of observers. The proper time integral of 
the force four-vector will be the momentum four-vector. We will not, here, 
discuss this formalism itself. Rather, we verify that this procedure gives 
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the Ehlers-Kundt result for plane-fronted gravitational waves, as may be ex- 
pected on account of the similarity of the two ideas. When we apply it to 
cylindrical gravitational waves it is found that the result so obtained is phys- 
ically reasonable and gives an exact expression for the momentum imparted 
to the rest particle, corresponding to the approximation given by Weber and 
Wheeler. 

In the free fall rest-frame the extended ?/>N (c0N) force four-vector is 
given [5] by 

F = m[(kiAf) fi + g ik g jkfl g jl g Ufi /AA], F t = m{\n^g~^)^ (1) 

where A = (In y/— g) tQ and g = det(gij). Thus the momentum four-vector, 
Pix, is 

ViX = J F^dt. (2) 

2 Plane-Fronted Waves 

The metric for plane-fronted gravitational waves is [6,7] 

ds 2 = dt 2 - dx 2 - L 2 (t, x){exp[2/3(t, x)]dy 2 + exp[-2(3(t, x)}dz 2 }, (3) 

where L and (5 are arbitrary functions subject to the vacuum Einstein equa- 
tions, 

L, aa + L(3% = 0, a = 0, 1. (4) 
Since L and (5 are functions of u — t — x, eqs.(4) reduce to the single euqation 

L uu + L(3l = (5) 

From eqs.(l) and (5) 

F = -m(L + (3 2 L)/L = 0, F 1 = 0, (6) 

where a dot denotes differentiation with respect to t. Consequently the mo- 
mentum four- vector becomes = constant. Thus there is a constant energy 
and momentum. The constant, here, determines the "strength" of the wave. 
This exactly coincides with the Ehlers-Kundt method in which they demon- 
strate that the test particles acquire a constant momentum and hence a 
constant energy, from the plane gravitational wave. 
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3 Cylindrical Gravitational Waves 

The metric for cylinder gravitational waves is [2,3] 

ds 2 = exp[2(^-^)]dt 2 -exp[2(^-^)}dp 2 - p 2 exp(-2^)d<p 2 -exp{2^j)dz 2 , (7) 

where 7 and ip are arbitrary functions of the time and radial coordinates, t 
and p, subject to the vacuum Einstein equations, 

i/>" + -i// -xj> = 0, i = p(^ ,2 + ^ 2 ), 7 = 2p^'- (8) 
P 

where a dot denotes differentiation with respect to time and a prime differ- 
entiation with respect to p. 

The solution of eqs. (8) is given by [3] 

ip = AJ (up) cos(cjt) + BY (cup) sm(ut), 

7 = icjp{(A 2 J Jo - B 2 YqY^) cos{2ujt) - AB[{J Y^ + Y^) sin(2wt) 
- 2(J Y^-Y J[ ) )ujt}}, (9) 

where A and B are arbitrary constants corresponding to the strength of the 
cylindrical gravitational waves, J and Y are the Bessel function and the 
Neumann function of zero order respectively. Here a prime denotes differen- 
tiation with respect to up, u being the angular frequency. 
From eqs.(l) 

F = -m{u[AJ cos(cjt) + BY am(ut)] - 2pu[(A 2 J J - B 2 Y Yq) cos(2cut) 

- AB(J Y^ + Y J' Q ) sin(2cut)] + 2[AJ sin(cjt) 

- BY cos(ujt)] 2 /AJ sin(cut) - BYb cos(cut) 

+ 2cjp[AJ sin(cjt) - £y cos(u;£)][A/ cos(u;i) + BYq sin(cut)]}, (10) 
Fi = m{AJ cos(cjt) + firj sin(cjt) - - [(A 2 J 4 ~ # 2 W) 
+ Lup(A 2 J 4 - B 2 Y Yq)'] cos(2c^), 

- iAB[2(Jo^' + I0J0) " + Y o4)'] sin(2^) 

- ^AB[A{J Y^-Y G J t Q ) + 2up{J^-Y G J' Q )']ut}. (11) 
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The corresponding p and momentum, p 1; imparted to the test particle is 

Po = — m[lnAJo sin(c<jt) — BY cos(cut) | + In |1 — 2up[AJ' cos(cjt) 

„ T r i / mi/ A 2 Jn JL -\- B^YoYn . 
+ ^0smM)]|(l + ^ 2 + ^ 2 ) 

AB(AJ Y Y ' + Y 4 - Alu P J J>Y{ - J Y>) 



x tan 



ujp^AJ' 2 + BYq 2 )^J1 - Auj 2 p 2 {A 2 J' 2 + B 2 Yf) 
1 (1 + 2AujpJ' Q ) tsm(\ujt) - 2BujpY Q ' 



l-Auj 2 p 2 (A 2 J' 2 + B 2 Y n ) 



AB(J' Y - JqFJ) 

+ Ml^T5^) +/lM] ' (12) 

777- I. 

Pi = -{[AJ' Q sin(c^) - Sy ' cos(wt)] " -M 2 Jo Jo ~ B 2 Y Yq) 
uj 4 

+ pu(A 2 J J' - B 2 YqYq)'\ sin(2cut) - \aB[J Y q ' + Y J' Q ) 

+ pcu(J y ' + Y Q J' Q )'\ cos(2cut) - ABuj 2 t 2 [(J Y ' - F J ) 

- puj(J Y^-Y J>)] + f 2 (puj)}, (13) 

where f\ and /2 are arbitrary constants of integration. Weber and Wheeler 
[2,3] exclude solutions that contain the irregular Bessel function, Y (up) as 
not well defined at the origin. Taking the Weber- Wheeler solution, eqs.(12) 
and (13) reduce to 

Po = — m[ln |AJ sin(a;t) | + (1 + AK /ujpJ' ) In |1 — 2ujpAJ' Q cos{ujt)\ 

+ /iML (14) 

?77 ! 

pi = —{AJ' sm{LJt) - -A 2 [J J' + pu(J J' Q )']sm(2ujt) + f 2 (poj)}. (15) 
uj 4 

We see that the quantity pi given by eq.(15) can be made zero for the small 
and large p limits by choosing f 2 — 0. This expression coincides with the 
approximate momentum expression given by Weber and Wheeler in the large 
and small p limits. Hence this is a physically reasonable expression for the 
momentum imparted to test particles by cylindrical gravitational waves. The 
quantity p given by eq.(14) remains finite for small p and can also be made 
finite for large p by choosing fi = — ln(J ). However, there is a singularity 
at out = mr. This problem does not arise in the general expression given by 
eq.(12). However, in that case there appears a term linear in time which 
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creates interpretational problems. Also p and p\ become singular at p = 
if B ^ 0. 

4 Conclusion 

The Ehlers-Kundt method gives the physically reasonable result that plane 
gravitational waves impart a constant energy and momentum to test parti- 
cles in their path. However, it does not provide a simple formula that can 
be applied to other cases. Since the essential idea there is exactly embodied 
in the e^N-formalism [5] it is reasonable to expect that it would give the 
same result. We have seen that the e^N-formalism gives a physically ac- 
ceptable expression for the energy and momentum imparted to test particles 
by gravitational waves which coincides with the Ehlers-Kundt approach for 
the plane-fronted waves. The expression for momentum coincides with the 
Weber- Wheeler [2,3] approximate result for cylindrical waves for small and 
large values of p. In effect it provides a general formula for the momentum 
imparted to test particles in arbitrary spacetimes. 

The inclusion of Y in the solution increases the energy at the source and 
thus provides a singularity at p = 0. If we do not include the Fo term in 
the solution, po becomes infinite at various times. One would have normally 
taken p to be the energy of the test particle. However, that should be 
given by (m 2 + pl) 1 ^ 2 . As such the significance of p is not clear. As this 
is the quantity which has the problem of becoming singular, while p\ is well 
behaved we are ignoring the problem. However, that problem does need to 
be resolved and an interpretation for po be provided. 

There is a problem of defining energy in general relativity which is most 
severe when trying to understand the energy content of gravitational waves. 
As gravitational waves are solutions of the vacuum field equations the stress- 
energy tensor is zero. This is the problem regarding the reality of gravi- 
tational waves. On the one hand it seems that gravitational waves do not 
carry energy but on the other hand detectors are built to extract energy 
from them. If we linearize the Einstein field equations [3,6], the complete 
non-linear equations can be rewritten with the linear part on the left side 
and all the non-linear terms on the right side. The non-linear part can then 
be regarded as the stress-energy tensor for the linearized wave equation [6]. 
Also, a 3 + 1 split, such as that of Arnowitt, Deser and Misner [8], provides an 
expression for the energy density in gravitational waves. We are only provid- 
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ing an expression for the energy imparted to a test particle. This expression 
is based on an operational 3 + 1 split. Hence this procedure might be able to 
provide a handle to the problem of the energy content in gravitational waves 
in particular and the definition of energy in general. To be able to apply it 
successfully the problem of p would have to be resolved. 
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